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1. Introduction 

There is strong observational evidence for a large scale galactic magnetic field in our 
galaxy and in neighbouring galaxies . From the early polarisation of starlight measurements 

01 [|2| to the more recent observations of the faraday rotation of pulsar emissions or other 
extra-galactic sources, there is a strong experimental evidence of such field. The galactic 
magnetic field is also required to explain why cosmic ray composition makes it appear 
that they have been through 30 to 300 time as much matter as there is between us and 
the galactic center. Alfven and Fermi (1949) first discussed this and predicted a field of 
10//gauss to confine the cosmic rays in the galactic plane, this is very close to the current 

2 — 3/igauss value for the coherent component. Current theories of this magnetic field 
creation require a seed field which is then amplified by a dynamo mechanism. 

In addition, more recent observations show that there seem to be a magnetic field on 
cluster size. Any theory of large scale structure should predict some seed magnetic field 
uniform on a galactic scale, the texture mechanism Qdoes this in a very natural fashion. 

The outline of this paper is as follow. In section 2, we will describe how one generically 
generates magnetic field with a model of cosmological perturbations. In section 3, we will 
briefiy describe the model of large scale structure formation by textures, and quote some 
of the results we will be using. In section 4 we will apply the formalism of section 2 to 
the texture model of structure formation and calculate the magnetic field it generates. In 
section 5 we will discuss the results and the implications for texture and compare to other 
candidate model of structure formation. 

1.1. Notation 

When discussing physics happening around matter-radiation equality, it is natural to 
refer scales to that era, so we will use a scale factor normalised at matter-radiation equality: 
a(t), wherea(teg) = 1 . Sometimes we will use the more common scale factor defined as 
-^ = Tfi-^- Spatial distances denoted by x , will always be in co-moving coordinates 
in units of Mpc. Xf,q denotes co-moving coordinates with a value that they would have 
at Matter-Radiation equality, x without the eq subscript denotes the normal co-moving 
coordinates, referenced to the present time, k will denote ^ and will be in co-moving 
units of Mpc~^; here again k^q will be referenced to equality, and k without a subscript 
will be referenced to the present. 
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Time for a co-moving observer ,t, will be given in units of Is . Where 1 ls= 2-\/2cr* ~ 
2.79 X 10^^ seconds We will convert between a and t using 

t = ^^((l + a)i -3^1 + ^+2) (1.1) 

and it's numerical inverse. Variables with a subscript i ( i.e. ai , Ti ) refer to the value 
of that variable at the epoch when a given texture just collapsed. 

1.2. CDM Model 

In all following discussions standard CDM will refer to a Cosmology with a scale 
invariant Harrison- Zeldovich (HZ) spectrum. For both CDM and the Texture model we 
will use a model with critical density = 1, with a Dark Matter component O ~ 0.95 
and a baryonic fraction of 5% . For the standard CDM model, we will assume a Hubble 
constant of 50km/s/Mpc, while for the textures model we will assume a Hubble constant 
of Tbkm/s/Mpc . For our purposes, the exact composition of the model makes little 
difference as long as most of the matter is Dark Matter. 

Throughout this article we will use DM when speaking about the dark matter, while 
baryonic and leptonic matter and photons will be refered to as matter. 

2. Magnetic Field Generation 

When thinking about magnetic field generation, it is useful to separate the physics 
in 2 phases. First a seed field, Bgeed, is created before a galaxy collapses. In the second 
phase the magnetic field is first amplified by fiux conservation during collapse and then by 
a dynamo process. 

The seed field creation ceases after the protons and electrons present in the primordial 
plasma recombine. The fiux lines are then frozen in the interstellar gas due to the large 
conductivity of the plasma. 

For an introduction to the dynamo theory see ]§][§[ [0 • The dynamo amplification 
mechanism occurs in galactic objects. The fiux is amplified through processes that makes 
fiux lines buckle and form loops arcing above the galactic plane . These loops twist 180°, 
rejoin and fall back on the galactic plane, increasing the fiux line density. For example, a 
cloud of gas can rise due to buoyancy forces, it then expands due to the lower density of gas 
above the galactic plane. As it expands, the coriolis force makes the cloud rotate, twisting 



the lines of flux as they are frozen in the cloud. The cloud then cools and falls back onto 
the galactic plane, carrying the loop back in the plane. This process increases the random 
field, Brms- The ensemble average of these random fields has a growing coherent mode. It 
is a topic of some controversy as to whether this yields an exponential or a linear increase 
of the magnetic field fl^ Jill . 

To get the current value of 3//gauss,with exponential dynamo , seed field values of 
only about 10~^^ to 10~^^gauss are required at the time of galaxy formation. On the 
other hand a linear dynamo requires fields of at least 10~^^/xgauss which are very hard 
to generate without delicate adjustement of large-scale structure generating models. For 



such an example in the infiationary scenario see [^| 

The creation of a seed field of even 10~^^ gauss is a non-trivial test for a theory to 
pass. Generically seed field creation can happen in 4 possibles eras. The first one is if the 
magnetic field has a primordial origin, for example |1^. The next two other eras involve 



dynamical phenomenas. The last period happens once the structure formation is well 



advanced. It involves magnetic field expulsed from powerful beamed radio sources. [|13 

We will not deal here with theories postulating a primordial of the magnetic field or 
those postulating the existence of powerful radio sources at very high redshift. 

2.1. Harrison scenario 

The second era for magnetic field generation is between matter-radiation equality and 
recombination. If we assume that there was some primordial vorticity, such as could happen 
in chaotic Cosmology, we can generate magnetic fields using to the following mechanism 



due to Harrison [1^1 1 15 



During that era the electrons motion is tied to that of the photons though the large 
thompson cross-section, aj. For the protons ,the proton's thompson cross section, crj, is 
much less important than the hadronique neutron-proton cross section. As a result while 
the electrons move with the photon fiuid, the protons together with the neutrons can have 
a different behavior. 

Implicit in the derivation of the Harrison mechanism is the existence of primordial 
vortices. We can nevertheless define a local scale factor r within the framework of the 
Homogeneous world model. The eddy follows the expansion and pr^ and pr'^ remain con- 
stant for matter and radiation respectively. So the angular momenta of the matter, pcor^, 
and the angular momenta of the radiation, prCOr^^, are separately conserved. Therefore 
Up (X -^ and u^ ^ u^ oc - . This means that electrons and protons will have different 



angular velocities, leading to an angular electric field and therefore a poloi'dal magnetic 
field. The value of this field will get frozen in at recombination. 

The detailled kinematical analysis can be found in Appendix 1. The final result is 

2t77 

BrecomUnaUon = ^w(Zr) = [-7 AQIQ~'^'^ gaUSs{hs))v^{Zr){hs)~'^ (2.1) 

e 
In an expanding Universe, the initial vorticity (^ = V x t7 ) perturbations decay as 

C(a) ~ 4 (2-2) 

This mean initial vorticity is washed away unless the universe was very chaotic when 
it started. Limits on the amplitude of primordial density perturbations by the COBE 
satellite rule out this possibility unless the universe was re-ionized. 

2.2. Harrison-like scenario 

The above Harrison scenario fails to pass the experimental tests because it depends on 
primordial vorticity. We will ressucite the Harrison mechanism, but rely on dynamically 
generated vorticity. This will yield the third period when magnetic field can be generated. 
The vorticity in the matter will be generated when the DM collapses early and acquires 
non-spherical components through tidal torquing. The acoustic waves in the matter, that 
form as the super-horizon density modes enter the horizon, will then undergo oblique 
shocks because they are falling in a non-spherical rotating potential. This mechanism will 
thus create the required vorticity in the matter. 

In most theories of Large Scale Structure (LSS) formation the Jeans mass before 
decoupling is larger than the mass within the Horizon. This is due to the photon pressure 
which keeps density perturbations initially at rest from growing . 

But LSS theories with topological defects generate density perturbations and velocity 
until very late times, unlike standard LSS scenarios which have all their density perturba- 
tion as initial conditions. 

In the texture model the density perturbations are not fixed as initial condition. The 
density perturbation are dynamically generated because of the texture collapse. We will 
see later that the regions around the center of a collapsing texture get a velocity kick 
toward that center. The radial component can be approximated by 



v = ec(l-^) (2.3) 

where r is the real radius, ti is the coUapse time of the texture and e is the only parameter 
of the theory and is normalised to COBE. This radial kick allows a texture seeded velocity 
kick to grow density even though the matter contained within it is lower than the Jeans 
mass before decoupling. 

The velocity kick is calculated for an ideal spherically symmetric texture collapse. 
Numerical textures simulations show that most textures are not perfectly symmetric and 
that a small azimuthal velocity is also imparted to the matter surrounding the texture. 
Typically there are non-radial components to the velocity of approximatively 1 -2 % of 
( p.3|) . This is not to say that a vorticity is created, it is a perturbation from a purely radial 
velocity kick due to the assymetric nature of real textures. 

We will not be using the small rms azimutal velocity from the initial collapse for our 
final result. Instead the initial assymetries in the dark matter, inside the volume of the 
texture, will be caused by the previous collapse of the smaller textures that were contained 
in that volume. We will see in section 4 how those assymetries cause the DM halo to 
gain angular momentum through tidal torques from neighbouring DM halos. We will see 
how this rotating DM halo will then transfer angular momentum to the matter before 
the recombination era. This allows some magnetic field generation for late textures by a 
mechanism similar to the Harrison scenario. 



3. Large-scale structure texture seeds 

The theory of topological defects as source of large scale structure was first treated 
by Kibble []T^, later Turok [^showed that global texture could occur naturally in GUT 



theories and produce interesting density perturbations in the early universe. (See []T^ [p!8|for 
review) 

An SU(2) symmetry seems most natural in term of a particle physics model. With a 
complex higgs field we could have a global SU(2) symmetry . For computational purposes 
it is easier to deal with a SO (4) 4-component isovector with real fields, so let us consider 
a GUT theory with a global S0(4) symmetry . 

/: = i(a^$)(a^$)-^($2-$^)2 (3.1) 




with 

1 / fTi„ 

(3.2) 

where the lagrangian is invariant under the global symmetry transformation $ -^ exp la^ 
with a a constant. In the ground state , $ develops an expectation value, < 0|$|0 > = 
($07 0, 0, 0), which breaks the symmetry of the lagrangian. . The 1-direction is arbitrary 
so that, early in the universe, fields not in causal contact will point in different directions. 
As their event horizon grows and they come in causal contact they may form a field 
configuration known as texture. This large size field configuration is unstable by derrick's 



theorem. |jT^ ^^ which states that no static soliton solution exists in 3 or more dimensions. 
So the energy stored in the gradients of the field configuration collapses to the center. 
At that point the energy density is such as to allow the texture field configuration to 
topogically unwind over the potential. This leaves the field aligned within the current 
horizon volume. As the horizon grows, other texture field configurations may appear and 
collapse. Each time a texture collapses, the energy of the field configuration produces a 
gravitational attraction toward the center of the texture. The resulting radial velocity kick 
can be described as 

v = ec(l-^) (3.3) 

where e = SttG^^ is normalised to COBE. There is a resulting azimutal velocity kick of at 
most 1-2% of this. The resulting velocity field does not have a net vorticity. This aximuthal 
component is simply due to the field configuration not being uniformly spherical. 

3.1. Power spectrum 



In [^ Pen et al.. calculated numerically the power spectrum caused by texture as it 
would look today. They give a fitting function 

VPMk) ^ 890iMpc^)k 

(27r)3 (1 + 14. xfc + llx A;i-5)2 ^'^- ^ 

We can use the fitted function of the numerical results in [|21[and evolve it backwards 
in time using linear theory and an time dependant cut-off for small k 

D{af 
Diar 



V ryKQQYfioyj Cl) i' -T (^ftcomouingj ^ ti/' „ \2 \ ) 



where kcomov < ^ and D(a) = (1 -f ^) is the growing mode solution of a density 

perturbation valid in both the matter and radiation era. |2 



4. Magnetic field by texture 

As we saw in section 2, we have to separate the problem of magnetic field generation in 
2 phases. In the first phase the texture provided the initial aspherical density perturbation. 
In the second phase angular momemtum will be provided by torquing by gravitational 
forces. Note that this process can be also applied to standard CDM. 

4.1. Angular Momentum 

Once the Texture collapse has seeded a velocity kick in the DM particles, our mecha- 
nism still require that the resulting particle distribution acquires a non-spherical distribu- 
tion. This rotation of the non-symmetric gravitational potential will be acquired through 
the tidal torquing picture. While in standard CDM the Dark matter cannot collapse this 
early, the tidal torquing picture of Peebles [^can still be applied. In what follows we 
will closely follow |J2^who computed the tidal torquing for a standard CDM model with 



gausssian spectrum. 

Consider a perturbation which is growing, we can look at the evolution of shells. A 
given shell will contain a fixed amount of matter and shells will not cross at least until 
they are collapsing faster than the Hubble expansion. The gravitational torque induced 
on a shell, by all the surrounding matter is. 

T{x) = -Gf p{s)sxV^is)d^s (4.1) 

J xEshell 

where 

p{s) = p,{l + 6{s)){l + e{s)) (4.2) 

All the non-sphericity of p is in e(s) for s much larger than the size of the texture (i.e. 
the other earlier textures that collapsed within this volume provide the asphericity ). 
$ is the gravitational potential, which can be expanded in spherical harmonics. 

^(^1 = E ^ E «^-(«)^^-(^' ^>' (4-3) 

1=0 m=-l 

where 

/"OO 

aim{x) = p, / Yi^{e,(t)){l + 6{s))e{s)s-'-H''s (4.4) 

J X 
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plugging (^4. 3D into (^4.1j)and using the fact that 

xxV{x'Yi^{e,(j))) =ix'LYim{e,<f>) (4.5) 

where 

L = -i{x X V) = x(L+ + L_) + iy{L- - L+)zL^ (4.6) 

is the angular momentum operator familiar from quantum mechanics, we obtain for 
the z component 

°° x' ^ f 

T^ = -iGMsheii X] 2/ 4- 1 S rnaim{x) I e{x)Yimdn (4.7) 

;=o m=-l -^ 

Let us introduce the multipole moments of the shell 

qimix) = f YC^{e,<P)s'pis)dh = ^Msheii JYC^iOAXm^ (4.8) 

Jshell 47r J 



where the final form is reached after substituting (^]^)and being left with the non- 
symmetric part. We can now rewrite r^ as 



r^ = -z47rG ^ — -y Y^ maimQL (4-9) 

/=0 m=-l 

If we retain only the quadrupole term {I = 2) we can compute the rms torque on a 



shell. 

(I^l) = 3(|r,|) 



.22 M 1 Dl 

= 3(— G')2 Y^ Y rnn{a2m{x)a*2^{x)q^^{x)q2n{x)) 
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with 



{a2mix)a2n{x)q2mi^)q2n{x)) = 

^'^hf' r^r^^l 1 1 1 ^2*,n(l)>^2m(2)y2m(3)y2;(eie2e3e4)ciOici02ci03ci04 

(4.11) 
the final form is reached under the assumption that the 4 point function can be written 



as 



(eie2e3e4) = (eie2)(e3e4) + (£163) + (e2e4) + (eie4)(e2e3)) (4.12) 



which is certainly true for a random gaussian process [22], and a reasonnable approx- 
imation in our case. Then after using the definition of the power spectrum in terms of the 
2 point function 




Yi^{l)Y2n{2){eie2)dQidQ2 



{2nf 



-ik-xi -y* 



d'kP{k) e-^'=-^^y;^(i)rfOi e-''-^-Y;^{2)dn2 



-ik- 



(4.13) 



,the expansion of a plane wave in sperical harmonics 



-ik-x 



A'KY.i'jiikx) J2 Y;^{e,<j>)Y2m{0k,phik) 



(4.14) 



1=0 



■m= — l 



and orthogonality relations one gets 
r{x) =< \T{x)f >2: 

where 



< almi^) >< qimi^)- < a2m{x)q2m{x) >' 



4 r 

< (llmix) >= J^^^shJ k^dkPik,T){j2{kx)f 



< aU^) >= ^^^^^^^VM!, I dkP{k,r){n{kx)f 



TT 



< a2m{x)q2m{x) >= — ^^"^ — — / kdkP{k,T)ii{kx)J2{kx) 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



we find that the volume factor exactly cancels out the volume factor of the power 
spectrum, for the sake of brievety, we omitted the r, r^ dependance. Putting it back, and 
integrating we find 



J(x, a, tti 



3t 



w{x, a, tti] 



eq 



da'a'ri^a' , xq, a^) 



M^heiix'^ ^v^,^,^v 2{2-^/2) J a^^,^,f VTT^x{a,ai)^M shell 
The M shell factor cancels exactly the Mass factor in (|4.15|) . f is given by 

f(xo,a',ai) r(x(xo,ai,a'),a') 



M 



M 



(4.19) 



(4.20) 
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and r is given by ( 4.15 ). x{xo, Ti, a') is the comoving position of a shell that was at 
position xq when the texture collapsed at conformal time ai . 

We can therefore use the texture power spectrum ( |3.5|) to figure out the angular mo- 
mentum produced by gravitational torquing. We also need the position of a shell as a 
functions of time. In the Texture model, Texture collapse before recombination, even in 
the radiation era. Collapse is modeled in the radiation era by 

x{xo,a,ai)=[ (4.21) 

V pb{a) I 



We have used the 5(a, ai) derived by |T^. When the shell enters the matter era, we match 



to the spherical collapse model [^by matching the peculiar velocity of the collapsed shell. 
Our integrals for the torque then runs until the turn-around time, as given by the spherical 
collapse model. 

The angular momentum thus calculated is not directly equivalent to vorticity. It 
can generate a solid body-like rotation of the dark matter fiuid spheroid but it has no 
net vorticity. To convert this rotation into vorticity, the fiuid should undergo dissipative 
processes or shocks. While the DM cannot undergo dissipative processes, it is still possible 
to convert the angular momentum from the DM shell into vorticity for the matter. 

As new modes of the density fiuctuations enter the horizon, the matter mode grows un- 
til it is pressure supported, and undergoes a few cycles of dissipatives acoustic oscillations. 
||22|| It is those shock waves, oscillating in the rotating aspherical gravitational potential 
that transfer the angular momentum to the matter. Obtaining the efficiency of this highly 
non-linear phenomena at transfering vorticity will await our future developpement of a 
non-linear hydro code to test the details of this model. Nevertheless we will make useful 
approximations. 



The torque is caused by a quadrupole-quadrupole interaction, whose strenght scales 



as ~ a ^, and whose amplitude, proportional to P(k), goes as ~ (1 -(- |a)^. The net result 



is that the torque is much more effective at early times. 

This is reffected in our results. Fig 1. shows the angular velocity acquired by an 
object of mass 2 x lO^^M© seeded at scale factor a^ = 0.032 at a comoving radius of IMpc. 
Most of the angular velocity of the DM is acquired early on. 
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angular velocity for a galaxy sized seed 



a/a_eq 



Fig. 1: angular velocity at a=50 for a texture collapsing at 0.032aeg the units 
are (Is^^) 1 Is = 2.79E12 seconds. The coinoving radius is 1.014 Mpc and 
this shell contains 2 x 10^^ Mq 

In Fig 2. we have plotted the angular velocity attained by each shell at their own turn 
around time. Note that the steep distribution at low radius is made steeper because by 
the time the outer shells have collapsed, the center will have collapse even further, causing 
a spin up if angular momentum is conserved. This could possibly slow the collapse of the 
center. 

In table 1 we have tabulated data. The angular velocity as a function of initial 
comoving radius follows roughly w ^ ^ and w ~ ^ we can compute the magnetic field 



for different size textures from ( |2.1| ). After the radiation era, w ~ aao so we can take w 
approximatively constant. As noted in the appendix we need the angular velocity of those 
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w at turn around time of each shell ai=0.032 



D.4 0.£ 



xO (Mpc) 



Fig. 2: angular velocity at the turn-around time as a function of the unper- 
turbed comoving radius for a texture collapsing at 0.032aeg. Note that each 
xta is reached at a different time, the units are (Is^^) 1 bs = 2.79 x lO^'^ 
seconds. A comoving radius of contains 2 x W^^ Mq(y^—^)^ 

shells that reached their turn around point before reconibination( a ~ 5). This gives that 
for a typical galaxy-sized object (aj = 0.048,xo = lMpc,Zf = 116, =^ Mgai = lO^^M© ) 
only regions inside a comoving radius of 0.2 Mpc have reached their turn-around radius. 
For those radii we obtain uj ~ 0.3/s~^. To compute a magnetic field we need to know 
how efficiently has the DM transfered it's angular momentum to the matter component. 
Let us get an upper limit to the possible magnetic field. Assuming 100% conversion, the 
matter acquires an angular velocity of w. We use (|2.1|) to get B ~ 2 x 10~^^ gauss. At 
that ealy epoch, the density of the proto galaxy is actually bigger than that of the final 
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object. Upon complete collapse to the density of a galaxy, this field will be multiplied by 
/_Pconaps^\2/^ ~ 3 X 10~^ So textures can yield at most a seed field of 6 x 10~^^ gauss for 
galaxy sized perturbations. 

Our galaxy has a rotational velocity of tw ~ 0.0005 in our units. The generated 
angular velocity for a galaxy-sized object was uj ~ 0.3/s~^, at recombination, which we 
have to multiply by the windup factor, 3 x 10~^, to yield 10"'^. This number has to be 
multiplied by the efficiency factor of conversion of DM angular momentum into matter 
vorticity. This shows that the required conversion factor cannot be as high as 100% .but 
has to be less about 10% . Under that assumption, the angular momentum of galaxy-sized 
objects produced by the texture scenario is adequate. 



5. Conclusion 

In the early universe, the texture model can produces sizable angular momentum. A 



similar analysys using the power spectrum of CDM has been peeformed by Ryden ||23|| . 
For both models the values are consistent with our own galaxy. 

In the texture model the early formation of angular momentum yields to formation 
of a magnetic field. On the other hand for the CDM model no magnetic field can be 
generated by that angular momentum because all of it was acquired after recombination. 
Our Scenario can only produce magnetic field in Clusters if a hiearchical model of cluster 
formation is adopted. This model only produced magnetic field in Galaxy sized objects 
and none in the interstellar medium other than that carried by the gas ejected from proto- 
objects. 

The upper limit we compute for the magnetic field falls short of the modern estimates 
of the required seed field. Although the texture model fares better in this problem than 
the standard CDM model, it does not pass this test. Magnetic field creation through early 
AGN is a better possibility in the texture model than in CDM because the structures form 
earlier. 

The scenario we have developped shows us that any model that yields non-baryonic 
growing density perturbations by the epoch of recombination can produce magnetic fields. 
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Table 1 








tti 


Xq 


<^/ Ojeq 


w (/s-i) 


0.003 


0.0139606 


0.134108 


20.6876 


0.003 


0.0698029 


1.11764 


3.46906 


0.003 


0.139606 


6 


1.25885 


0.004 


0.0184639 


0.18984 


15.9377 


0.004 


0.0923195 


1.58212 


2.33627 


0.004 


0.184639 


8 


0.773242 


0.006 


0.0274655 


0.31313 


10.6319 


0.006 


0.137327 


2.60968 


1.28811 


0.006 


0.274655 


12 


0.334348 


0.008 


0.0364617 


0.450027 


7.92929 


0.008 


0.182309 


3.75069 


0.789084 


0.008 


0.364617 


16 


0.167428 


0.012 


0.0546771 


0.762122 


4.83116 


0.012 


0.273386 


6.35207 


0.33716 


0.012 


0.546771 


24 


0.0595082 


0.016 


0.0728387 


1.1155 


3.26369 


0.016 


0.364193 


9.2977 


0.167547 


0.016 


0.728387 


32 


0.0288167 


0.024 


0.109059 


1.93168 


1.76581 


0.024 


0.545294 


0.805078 


0.057496 


0.024 


1.09059 


48 


0.0101015 


0.032 


0.144989 


2.87665 


0.912561 


0.032 


0.724945 


1.19898 


0.0274236 


0.032 


1.44989 


64 


0.00471943 


0.048 


0.216372 


5.12164 


0.287895 


0.048 


1.08186 


42.6972 


0.0102267 


0.048 


2.16372 


96 


0.00157769 


lis table we lis 
shells turns a 


t the angular ve 
round. This va 


ocity acquired 
ue is calculatec 


by the DM shell i 
i at the epoch o 



Int 
the DM 

each shell. The scale factor at that turn around time is given in column 3. These quantities 
are given for a range of initial texture collapse and initial co-moving radius. 
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Table 1 








tti 


Xq 


S'/Q-eg 


w (/s-i) 


0.064 


0.28702 


7.80306 


0.129639 


0.064 


1.4351 


65.0556 


0.00479661 


0.064 


2.8702 


128 


0.000734812 


0.096 


0.426246 


14.3865 


0.0423705 


0.096 


2.13123 


119.955 


0.00165626 


0.096 


4.26246 


192 


0.000283749 


0.128 


0.562777 


22.5143 


0.0190206 


0.128 


2.81388 


187.74 


0.00081579 


0.128 


5.62777 


256 


0.00015816 


0.172 


0.746374 


36.1093 


0.00949166 


0.172 


3.73187 


301.128 


0.000447707 


0.172 


7.46374 


344 


9.5687e-05 


0.256 


1.08455 


69.5543 


0.00365715 


0.256 


5.42276 


512 


0.000228605 


0.256 


10.8455 


512 


5.06313e-05 


0.384 


1.57217 


138.656 


0.0014671 


0.384 


7.86085 


768 


9.83743e-05 


0.384 


715.7217 


768 


2.08143e-05 


0.512 


2.03084 


228.674 


0.000952099 


0.512 


10.1542 


1024 


5.26307e-05 


0.512 


20.3084 


1024 


1.04759e-05 


0.768 


2.87671 


467.535 


0.000655441 


0.768 


14.3835 


1536 


2.13841e-05 


0.768 


28.7671 


1536 


3.3304e-06 


1.024 


3.64583 


779.061 


0.000399587 


1.024 


18.2291 


2048 


1.0869e-05 


1.024 


36.4583 


2048 


3.00253e-07 


1.536 


5.01273 


1598.52 


0.000185055 


1.536 


25.0637 


3072 


3.60041e-06 


1.536 


50.1273 


3072 





2.048 


6.21222 


2655.68 


0.000103789 


2.048 


31.0611 


4096 


9.4301e-07 


2.048 


62.1222 


4096 





2.99 


8.12198 


5168.26 


4.47294e-05 


2.99 


40.6099 


5980 





2.99 


81.2198 


5980 
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7. Appendix 1 

Here we will derive the equation for the Harrison mechanism. 

is the standard convective derivative where ^ follows the expansion and Vj = Rx + Rx. 
Rx is the hubble flow with the other term being the the peculiar velocity. In an expanding 
Universe it must be replaced by 

where we often replace ^= H , which is the hubble expansion rate. 
The equation for the electron motion is 

^;^(^^e) = —[E+-V,^B--) + -p^n,caTAv^-v,)- — -^- (7.3) 

the ion equation is 

f ;^(^^p) = +^(^ + c^^ "" a^ + -p-,n,caT,{v^ _ ^^) _ _J1 _ ^ (7.4) 

and the photon equation is 



1 ^1 D ,^. 1 -. rf , 1 ,„, 4 ,^ ^, 4 /^ ^^ 

(p^+^p^)- — (i^t'^j = --^V^—p^-{p^+-^p^)V(t)--p^neCaTS'"^-Ve)--p^npCaT^{v^-Vp] 

(7.5) 

+ \^^'v, (7.6) 

6 bna 

We drop the last term, because the conductivity is very large at that epoch. If we use the 

continuity equation 

V.V= , -\ , ^ (7.7) 



(Pt + ^Pi) dt 



and p^ = kp-yC^ and V ■ Vj = 3H + V ■ v 



3' 

we get 



P-y-p^v^ = -T^P-yy-yC^ ■ Vj) - PjVcj) - car^p^n^iv^ - v^) - car^p-yUpiv^ - Vp) {7.i 
JJt 6 
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Note that we will be taking the curl of ([f73| ) ( |7.4|) and (|77^)so that the pressure and 
gravitational term will go away. A short calculation of Vx( [7.4|) yields 



' ^ iR^C)+C{^-v)HC-^)v = -^^^{R'B) + -^Uv'B-^'B^ 



R^ Dt PpC R"^ Dt 47r(Jcond I /i c'^ I 3 Pp 

(7.9) 
where to find the RHS we assumed that B ■ W = HB which is only exactly true at 
9 = ^ . The conductivity a being very high, we can drop the wave equation term in B . 
With 



n„e 



«p = ^ = — (7.10) 

PiC rUpC 



we can finally write V x (|7.4| ) as 



1 D 



similarly Vx ([7. 3|) yields 

1 D 



[R^{Cp + apB)]=-\P:LcarAC,'Qe) (7.11) 



[R\Ce - a J)] = ^^I^^(C; - Ce) (7.12) 



RDt^ ^ •= /J 3 ^ 



e 



Right away we notice that the electron dynamics is much more affected by the photon 
drag due to the smaller electron mass on the denominator of the RHS of ( [7.12 ). Also 



notice that in ( [7.11|) and ( [7.12|) cry^/mp = ar /riip. Taking the curl of (|7.8|) yields 



P^ — — {RCj) = -neCaT,pgamma[K-i " Ce " Cp) (7.13) 



where we have assumed V ■ iT is small . 

((E1D+ (E33)) yields 



\^,m^. + C,)]=~[R'(c,B)] (7.14) 

We have used a^^ cxp- Since the proton will spin down faster than the electron, the 
proton vorticity will be small. Integrating ( [7. 14] ) yields the wanted result ( |2.1| ). 

This was derived under the assumption that the matter background is spinning down 
faster from having the same initial angular velocity as the radiation. In the original 
Harrison mechanism, the vortices were primordial and the radiation had the same initial 
velocity as the matter. The evolution was different because for the radiation (and the 
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electrons tied to it by the strong Thompson cross section) L = pVr uj ~ a~ a a uj = 
constant so uj-^ ~ a~^ where a is the scale factor, whereas for the matter w ~ a~^ . 

In our model the angular velocity will be provided by angular torquing (+ dissipation). 
We have to require that the angular momentum gained by the matter has had time to spin 
down. For angular velocity acquired by gravitational torquing ( + dissipation) , most of 
the angular momentum of the matter should be acquired by the time the dark matter 
shell has reached it's turn around radius. Before recombination the photon mean free path 
l-y = ^ ^^ is 6 X 10~^ MpcflQ^ h'^Q (ytj)^ so that density perturbation of the radiation by 
textures don't have time to diffuse out of their creation region before recombination and 
they will acquire angular velocity to due gravitational torquing+dissipation just as the 
matter did. So at a given shell radius, the angular momentum growth must peak before 
recombination ( a ~ 5 ) in order to have a sizable magnetic by the Harrison mechanism. 
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